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Abstract 

Quantum-chemistry  calculations  were  carried  out  to  predict  not  only  the  structure,  but  also  the  energetic  and  electrical  parameters  of 
lithium  trivanadates  of  the  formulae  Lii+AQOs  (where  x  =  1,2  and  3).  Altogether  21  entities  were  considered:  four  with  the  formula  LiV3Og, 
14  with  the  formula  Li2V30g  two  with  the  formula  Li3V30g  and  one  with  the  formula  Li4V30g.  A  lithium-vanadium  bronze  with  the  formula 
Li,  3V30g  was  synthesised  and  investigated  in  order  to  compare  experimental  and  predicted  characteristics.  The  presence  of  Li]+xV3Og  units 
in  the  synthesised  material  was  confirmed  by  infrared  (IR)  measurements.  The  predicted  value  of  the  open  circuit  voltage  (OCV)  change  for 
the  insertion  of  one  lithium  atom  equal  to  0.559  V  compares  well  with  the  experimental  value  of  0.52  V  for  the  insertion  of  the  first  lithium 
atom  into  the  synthesised  bronze.  These  results  form  the  basis  for  predicting  the  properties  and  modelling  or  technologically  interesting 
electroactive  materials.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

More  than  20  years  have  passed  since  the  lithium-vana¬ 
dium  bronze  (Lii+^VgOg)  was  used  as  an  alternative  to 
V6013  cathode  material  in  secondary  lithium  cells  [1], 
The  crystalline  Lix .iV3Og  is  monoclinic  and  belongs  to 
the  P2\lm  space  group.  It  is  composed  of  two  basic  structural 
units,  namely,  a  V06  octahedron  and  a  V05  distorted 
trigonal  bipyramids  [2]. 

It  was  found  that  the  electrical  properties  of  Li1+jrV308, 
such  as  discharge  capacity,  discharge  rate  and  reversibility 
and  cyclability,  depend  on  the  method  of  synthesis  of  the 
compound  and  the  conditions  under  which  the  cathode  is 
prepared.  Several  methods  of  production  and  processing 
Li  |  +  ,-V3Ox  electrodes,  have  been  reported  they  are  based 
on:  fast  cooling  [3],  hydrothermal  treatment  [4],  substitution 
of  lithium  with  other  monovalent  cations,  (like  Na+  or  K+) 
[5],  dehydration  of  the  lithium  vanadate  aqueous  gel  [6], 
intercalation  of  inorganic  molecules  (NH3,  H20,  C02)  into 
the  interlayer  [7],  ultrasonic  treatment  [8], 

The  electrochemical  properties  and  insertional  abilities  of 
Li1+YV3Og  cathode  depend  on  their  structure.  It  is  known 
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that  amorphous  bronzes  have  higher  discharge  capacity  than 
crystalline  ones  [9].  Kawakita  et  al.  have  confirmed  earlier 
reports  that  the  insertion  rate  of  the  lithium  ions  depends 
strongly  on  the  degree  of  crystallinity  of  the  material  [10]. 
The  electrochemical  and  physico-chemical  properties  of  this 
compound  have  also  been  studied  by  other  researchers  using 
EPR,  7Li  nuclear  magnetic  resonance,  as  well  as  infrared 
(IR)  and  Raman  spectroscopies  [11-13], 

Application  of  quantum  and  statistical  mechanics 
becomes  an  useful  tool  in  predicting  the  structure  and 
electrochemical  properties  of  electroactive  materials.  So, 
the  value  of  open  circuit  voltage  (OCV)  changes  for  the 
lithium-MnOo  (spinel  or  globular)  cell  was  predicted  at  the 
ab  initio  level  using  a  pseudopotential  [14].  With  method, 
the  Gibbs  free  energy  changes  accompanying  Li  |  _,Co02 
formation  were  calculated  for  different  number  of  Li+-ions 
in  the  host  structure  [15].  Computer  simulation  methods 
were  used  to  analyse  the  initial  stages  of  Li+  migration 
across  the  V205/Li2Si03  phase  boundary  [16],  The  local- 
density  functional  theory  was  also  used  to  study  the  struc¬ 
ture  and  phase  transitions  of  the  monoclinic  Li1+YV308 
[17]. 

In  this  work,  we  show  how  the  ZINDO  quantum  che¬ 
mistry  method  can  be  used  to  model  the  structure  of 
lithium  trivanadate  and  predict  its  electrical  and  energetic 
parameters.  We  further  demonstrate  how  the  predicted 


0378-7753/02/$  -  see  front  matter  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
PII:  S0378-7753(01)00986-7 


62 


E.P  Koval’ chuk  et  al.  / Journal  of  Power  Sources  107  (2002)  61-66 


chracteristics  compare  with  the  experimentally  determined 
structural  (IR)  and  electrical  (coulometric  titration)  charac¬ 
teristics  of  this  electroactive  material. 


2.  Experimental 

The  lithium-vanadium  bronze  was  prepared  hydrother- 
mally  from  LiOH  and  V205  at  a  temperature  of  about  50  °C, 
according  to  the  method  described  in  the  literature  [9].  The 
synthesised  bronze  was  washed  several  times,  first  with 
distilled  water,  then  with  methanol,  after  which  it  was 
vacuum  dried  at  150  °C  for  24  h.  The  cathode  was  a  mixture 
of  the  lithium  trivanadate  thus  obtained,  an  electroconduct¬ 
ing  additive,  and  PTFE  powder  in  the  weight  ratio  3:2:1, 
ABG-74  expanded  graphite  (Superior  Grafite  Company, 
USA)  or  TA-80  acetylene  soot  (USSR)  served  as  electro¬ 
conducting  additives.  After  the  addition  of  acetone,  the 
mixture  of  the  above  three  components  was  kept  in  suspen¬ 
sion  for  2  h  and  subjected  to  ultrasonic  treatment  for 
1 .5  min.  After  evaporation  of  the  acetone,  the  mixture 
was  vacuum  dried  at  200  °C  for  24  h  and  subsequently 
pressed,  at  a  pressure  of  5  MPa,  into  a  nickel  net.  Prismatic 
cells  were  then  produced  in  which  the  cathode  was  the  one 
described  above  and  the  anode  was  lithium  foil. 

The  content  of  the  lithium  ions  intercalated  into 
Li  i^V^Ox  was  determined  by  coulometric  titration  at  a 
current  of  300  pA  applied  for  300  s  as  described  before  [18]. 
IR  spectra  of  Li|+AV308  dispersed  in  KBr  (4000-400  cm  1 
region)  or  polyethylene  (500-100  cm-1),  were  recorded 
using  a  Bruker  IFS  66  FTIR  instrument  with  resolution 
varying  between  1-4  cm-1. 

3.  Computational  procedure 

The  structure  and  physico-chemical  characteristics  of 
the  lithium  trivanadates  were  predicted  at  the  level  of  the 
ZINDO  semi-empirical  quantum  chemical  method  [19] 
included  in  the  HYPERCHEM  program  package  [20]. 
The  21  most  probable  and  energetically  favourable  entities 
originating  from  Li1+YV308  were  initially  optimised  using 
molecular  mechanics.  These  structures  were  then  subjected 
to  full  quantum-mechanical  calculations  in  order  to  obtain 
final  characteristics. 


4.  Results  and  discussion 

Fig.  1  depicts  four  different  entities  of  formula  Li  V308 
containing  pentavalent  vanadium,  which  represent  the 
classical  views  of  the  bonding  of  elements.  From  the 
chemical  point  of  view,  the  lithiation  of  Li  V308  involves 
the  electrochemical  reduction  of  vanadium:  the  insertion  of 
the  first  lithium  atom,  leads  to  the  reduction  of  one  vana¬ 
dium  atom  from  +5  to  +4,  insertion  of  the  second  lithium 
atom  leads  to  the  reduction  of  the  next  vanadium  atom;  and 
so  on.  Pentavalent  vanadium  exists  in  the  singlet  state, 
whilst  the  quadrivalent  vanadium  is  in  the  doublet  state. 
Therefore,  during  the  intercalation  of  one,  two,  and  three 
lithium  atoms  the  multiplicity  of  the  lithiated  bronze 
changes  to  doublet,  triplet  and  quartet,  respectively.  In 
consequence,  the  number  of  possible  forms  increases  sub¬ 
stantially:  at  least  14  entities  of  formula  Li2V308  should  be 
considered  (Fig.  2). 

The  initial  stage  of  lithium  insertion  into  LiV308  was 
investigated  experimentally  and  the  analysis  indicated 
which  Li2V308  forms  are  most  favourable.  In  the  case  of 
higher  lithiated  bronzes  Li , +vV3Ox  (where  x  =  2  and  3), 
only  three  forms  were  taken  into  consideration  (Fig.  3). 

The  bond  energy  of  a  compound  (£]„„)  can  be  considered 
a  criterion  of  its  stability.  (Lhin)  is  the  difference  between  the 
energy  of  the  compound  and  the  energy  of  two  atoms  at 
absolute  zero.  According  to  the  molecular  orbital  theory,  the 
bonding  of  two  atoms  takes  place  when  electrons  move  from 
the  atomic  orbitals  to  the  bonding  molecular  orbital  which  is 
accompanied  by  an  decrease  in  energy  [21].  If  this  energy 
decrease  does  not  occur,  a  molecule  is  not  formed,  as  in  the 
case  of  2b  and  4c.  On  the  other  hand,  the  most  stable  entity  is 
3b,  for  which  Ebin  =  —2182.1  KJ/mol.  Other  possible  forms 
of  Li2V208  can  be  arranged  as  shown  below,  in  decreasing 
order  of  bond  energy: 

3b  <  2c  <  3c  <  3d  <  2d  <  2a  <  3a  <  lc  <  4a  <  4b  <  la 

(1) 

Doubly  lithiated  bronze  can  exist  in  the  two  forms  5a  and 
5b  presented  in  Fig.  3,  the  multiplicity  of  which  is  3. 
However,  according  to  quantum  chemical  calculations,  only 
5b  is  stable. 

In  the  case  of  trilithiated  bronze  the  vanadium  atoms  can 
all  be  quadrivalent  only  in  structure  6.  Therefore,  the 


Fig.  1.  Possible  structures  formulae  of  LiV308. 
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Fig.  2.  The  structural  formulae  of  the  monolithiated  bronze  (Li2V308). 


lithiation  process  can  be  presented  as  a  transition  between 
the  following  forms  of  the  bronze 
1  ->  la  ->  5b  ->  6  (2) 

The  energetic  characteristics  and  relative  charges  on 
selected  atoms  for  some  of  the  compounds  are  listed  in 
the  Table  1 .  The  increase  in  the  number  of  inserted  lithium 
atoms  causes  a  decrease  in  the  energy  of  the  lowest  unoc¬ 
cupied  molecular  orbital  (Z?lumo)-  This  indicates  that  with 
the  increase  in  the  degree  of  lithiation,  the  susceptibility  of 
the  bronze  to  reduction  increases.  On  the  other  hand,  there  is 
no  correlation  between  the  degree  of  lithiation  and  the  A Eq 
OElumo  —  /-homo )  gap,  °r  the  relative  charges  on  lithium  or 
the  relative  charges  on  lithium  or  vanadium  atoms. 


The  enthalpy  change  for  bronze  lithiation  is  roughly  equal 
to  the  Gibbs  free  energy  change  (AGr)  since  the  entropy 
change  for  this  process  is  insignificant  and  can  be  neglected 
[14]. 

The  OCV  changes  for  lithium  insertion  (A E)  can  be 
calculated  from  the  Gibbs  free  energy  changes  in  accordance 
with  the  equation. 


A  E  = 


-AGr 

(x2  —  X\)F 


(3) 


in  which  ( —  X\ )  denotes  the  stoichiometric  coefficient  for 
the  intercalation  of  Li+  into  the  bronze  cathode;  and  F  is  the 
Faraday’s  constant.  Table  2  sets  out  the  Gibbs  free  energy 
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Fig.  3.  The  structural  formulae  of  the  bilithiated  (Li3V308)  and  trilithiated  (Li4V308)  bronzes. 
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Table  1 


Energetic  characteristics  and  relative  charges  on  atoms  in  Li1+;tV308 


Compound 
(formulae  no.) 

^bin 

(kJ  x  mol-1) 

\Hf 

(kJ  X  mol~' 

^HOMO 

(eV) 

^LUMO 

(eV) 

A  Eq 
(eV) 

Relative  Mulliken  charges 

V(I) 

V(II) 

V(HI) 

Li(I) 

Li(II) 

Li(III) 

Li(IV) 

LiV3Os  (1) 

-345.8 

3353.5 

-7.355 

0.228 

7.583 

0.53 

0.66 

0.80 

0.39 

Li2V3Os  (la) 

-560.7 

3299.5 

-6.410 

0.195 

6.605 

0.69 

0.84 

0.49 

0.23 

0.40 

Li3V308  (5b) 

-666.5 

3352.9 

-9.881 

-0.760 

9.121 

0.82 

0.73 

1.04 

-0.42 

-0.12 

-0.07 

Li4V308  (6) 

-720.6 

3457.8 

-8.930 

-1.042 

7.888 

0.88 

0.50 

1.07 

-0.33 

0.30 

0.50 

-0.24 

Table  2 

Gibbs  free  changes  and  OCV  value  for  the  intercalation  of  Li+  ions 


Intercalation  process 

AG  (kJ  x  moF1) 

E  (V) 

LiV308  4-  Li+  +  e — 

>Li2V308 

-54.00 

0.559 

Li2V308  +  Li+  +  e- 

->Li3V308 

53.35 

-0.552 

Li3V308  +  Li+  +  e — »Li4V308 

104.85 

-1.086 

changes  and  the  OCV  changed  corresponding  to  the  inser¬ 
tion  of  one  two,  and  three  Li  atoms  as  calculated  by  quantum 
chemical  methods. 

Clearly,  the  value  of  A Gr  for  the  intercalations,  of  the 
first  Li+  into  the  bronze  cathode,  and  F  is  the  Faraday 
constant.  Table  2  sets  out  the  Gibbs  free  energy  changes 
and  the  OCV  and  changes  corresponding  to  the  insertion  of 
one,  two,  and  three  Li  atoms  as  calculated  by  quantum 
chemical  methods. 

Clearly,  the  value  of  AGr  for  the  intercalation  of  the  first 
Li+  is  negative,  which  means  that  the  process  is  thermo¬ 
dynamically  favourable  and  should  proceed  spontaneously. 
On  the  other  hand,  intercalation  of  the  second  and  third 
lithium  atoms  is  a  forced  process.  Intercalation  of  lithium 
into  the  synthesised  lithium-vanadium  bronze  (0  <  x  <  1 ) 
was  investigated  by  coulometric  titration  (Fig.  4).  The  OCV 
changes  versus  x  were  measured  for  electrodes  containing 
either  expanded  graphite  or  acetylene  soot  as  a  electrocon- 


Fig.  4.  The  OCV  changes  in  the  Li1+;cV3O8/1.0  N  U1CIO4  +  propylene 
carbonate/Li  cell  vs.  the  intercalated  lithium  stoichiometric  coefficient 
(at)  for  two  different  electroconducting  additives-acetylene  soot  (1)  and 
expended  graphite  (2). 


ducting  additives.  The  differences  in  OCV  corresponding  to 
x  =  1  are  equal  to  0.52  and  0.40  V,  respectively.  These 
values  correlate  reasonably  well  with  those  predicted  the¬ 
oretically  (Table  2)  [21],  The  value  of  the  OCV  changes 
reported  by  Kawakita  et  al.  for  the  oxygen-deficient  lithium- 
vanadium  bronze  (Li1+YV308_,5)  are  very  close  to  our  data. 

As  was  mentioned  above,  four  entities  with  the  formula 
Li  V3Ox  are  possible  (Fig.  1),  and  in  each  of  them,  vanadium 
atoms  are  quadrivalent  and  exist  in  the  singlet  ground  state. 
However,  these  formulae  are  essentially  different:  in  the  first 
and  second  entities,  a  six-membered  ring  appears  and  the 
number  of  V=0  groups  changes  from  one  (in  1)  to  four  (in 
3).  Both  V-O-V  and  V=0  fragments  were  actually  found  in 
the  synthesised  bronzes,  and  the  presence  of  three  different 
V=0  bonds  of  lengths  0.1586,  0.1584  and  0.1583  nm  was 
confirmed  by  X-ray  analysis  [2],  These  latter  bonds  in 
monolithiated  bronze  can  be  assigned  to  V(I11)=0  (in  1), 
and  V(II1)=0  and  V(I)=0  (in  4),  since  1  and  4  are  energe¬ 
tically  the  most  favourable.  The  above  facts  undoubtedly 
imply  that  the  entities  we  have  considered  are  fragments  of 
real  structures. 

The  IR  spectrum  of  the  synthesised  bronze  (Lit .3V3O8) 
exhibits  bands  at  1010,  994  and  972  cm-1,  which  can  be 
assigned  to  the  stretching  vibrations  of  the  three  above  types 
of  V=0  groups  (Fig.  5)  Similar  results  were  reported  by  Dai 
et  al.  for  bronze  synthesised  at  300-350  °C  [22],  The 
absorption  at  734  cm-1  can  be  related  to  the  V-O-V  stretch¬ 
ing  vibration.  In  the  case  of  the  oxygen-deficient  bronze 
(Li1+A.V308_^)  the  relevant  peak  was  observed  at  750  cm-1 
[18],  The  absorption  between  600-400  cm-1  most  probably 
arises  from  V-O-V  and  V=0  bending  vibrations,  while  that 
at  295,  212,  161  and  86  cm  1  can  be  related  to  the  deforma¬ 
tion  of  the  polynuclear  chain  consisting  of  Li-O-fragments 
[13],  Lastly,  the  absorption  in  the  region  3550-3350  cm-1 
and  that  at  1615  cm-1  are  probably  due  to  Li-O-stretching 
vibrations. 

Theoretical  calculations  (Table  3)  indicate  that  modes 
at  high  wave  numbers  (the  region  between  3550  and 
3350  cm”1)  correspond  to  Li-O-stretching  vibrations.  For 
bronzes  with  different  degrees  of  lithiation  the  intensity 
of  bands  between  3400-3330  cm”1  increases  somewhat 
when  the  number  of  intercalated  Li  atoms  increases.  This 
absorption  shifts  to  2996  cm”1  if  x  =  0,  When  x  =  3,  the 
band  splits  into  two  peaks.  An  increase  in  the  degree  of 
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Fig.  5.  The  IR  transmission  spectra  of  the  initial  lithium-vanadium  bronze 
l.i  [  ,V3Ox.  (r~0,3)  in  the  (a)  4000  500  cm  and  (b)  100  650  cm  1 
ranges. 


Table  3 

Theoretically  predicted  wave  numbers  (in  cm-1)  of  vibrational  modes  of 
Li,V308 


Li1V308  (1) 

Li2V3Os  (la) 

Li3V3Os  (5b) 

Li4V308  (6) 

2996  s 

3405  w 

3331  m 

3400  w 

1241  s 

1522  w 

1519  w 

3373  m 

831  m 

1123  m 

1 153  m 

1501  m 

781  w 

1023  s 

709  s 

1083  m 

624  m 

734  s 

596  w 

1047  s 

483  w 

610  m 

595  s 

886  s 

345  w 

506  s 

341  w 

873  w 

321  w 

474  s 

302  w 

824  m 

212  w 

422  s 

160  w 

597  s 

338  m 

132  m 

481  s 

283  m 

114  m 

450  m 

261  s 

88  m 

374  s 

178  s 

48  m 

302  m 

102  m 

234  m 

85  m 

212  w 

157  m 

107  m 

lithiation  enhances  low-frequency  modes  which,  as  men¬ 
tioned  above,  correspond  to  the  deformations  of  the  poly¬ 
nuclear  chain. 

Therefore,  the  proposed  forms  of  Li1+jrV3Og  are  a  good 
reflection  of  the  actual  interatomic  bonding  in  lithium- 
vanadium  bronzes.  It  should  be  noted,  however,  that  the 
real  bronzes  are  infinite  structures  made  up  of  fragments  of 
the  entities  we  have  been  considering. 


5.  Conclusions 

The  structure  and  properties  of  lithium  trivanadates  were 
examined  with  respect  to  the  model  molecules  constructed 
following  the  classical  valence  theory.  Four  entities  with  the 
formula  LiV3Og,  14  with  Li2V30g,  two  with  Li3V30g  and 
one  with  Li4V30g  were  selected  for  theoretical  examination. 
The  energetic  characteristics  of  some  of  these  entities  were 
used  to  estimate  OCV  changes  during  the  lithiation  of 
Li1+YV3Og.  For  the  insertion  of  the  first  lithium  atom,  the 
calculated  value  of  the  OCV  change  is  0.559  V,  while  the 
experimental  values  are  0.52  and  0.40  V  for  electrodes  filled 
with  expanded  graphite  and  acetylene  soot,  respectively.  The 
correlation  between  the  calculated  and  experimental  values 
of  the  OCV  changes  for  the  lithiation  process  was  also  good 
when  the  second  and  third  lithium  atoms  were  inserted  [10]. 
The  results  obtained  for  the  model  molecules,  form  a  con¬ 
venient  basis  for  modelling  the  bronze  clusters  (containing 
64  and  128  atoms).  This  is  one  of  the  issue  on  which  we  are 
currently  focusing  our  attention. 
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